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We show that the ground state of an attractive Bose gas in 
a double well evolves from a coherent state to a Schrodinger 
Cat like state as the tunneling barrier is decreased. The latter 
exhibits super- fragmentation as spin-1 Bose gas with antifer- 
romagnetic interaction, which is caused by the same physics. 
We also show that the fragmented condensates of attractive 
and repulsive Bose gases in double wells lead to very different 
interference patterns. 

Since the discovery of Bose-Einstein condensation 
(BEC) in atomic gases, there has been many studies of 
the phase coherence properties of Bose systems. A fre- 
quently studied example is a repulsive Bose gas in a dou- 
ble well, where the ground state changes from a coherent 
state to a Fock (or number) state as tunnelling is de- 
creased 0] Q |^]. In a coherent state, the condensates 
in the two wells are phase coherent; whereas they are 
incoherent in a Fock state. A number of authors || || 
H have shown that despite their differences, coherent 
states and Fock states can not be distinguished by inter- 
ference experiments which measure the density profile of 
two overlapping condensates that are initially separated. 
The reason turns out to be a subtle one: The projective 
nature of the measurement process introduces coherence 
into a Fock state, so much so that it gives rise to the 
same interference pattern as the coherent state. In con- 
trast, there are little studies of the phase coherence of 
Bose gases with attractive interactions. Such studies are 
not only important for conceptual reasons, but also for 
current experiments 0. 

We shall see, both repulsive and attractive interactions 
will lead to "fragmented" condensates, (meaning more 
than one macroscopic eigenvalue in the single particle 
density matrix). However, the ground state of an attrac- 
tive Bose gas in double wells is in general Schrodinger- 
Cat like and is "superfragmented" 0. The latter means 
that the internal number fluctuation AiV 2 is enormous, 
of order ~ TV 2 . In contrast, AN 2 ~ for the Fock state 
(which is fragmented), and is of order N for coherent 
states. Because of this huge fluctuation, the measure- 
ment process is unable to project the system into a single 
coherent state as in the repulsive case, but instead a sta- 
tistical mixture of Fock-like states with unequal number 
of particles in each well. Consequently, the interference 
pattern of attractive and repulsive Bose gases in double 
wells are very different in the fragmented regime. 

Before proceeding, we shall comment on two recent re- 
lated developments. There have been increasing efforts 
in recent years to create a Schrodinger cat state in con- 
densed matter and atomic systems. The Schrodinger Cat 
state (or "Cat state" for short) is a superposition of two 
macroscopic quantum states. It is often used to illus- 



trate the peculiarity of Quantum Theory, which admits 
such superpositions even though they have never been ob- 
served on the macroscopic scale. The usual explanation 
is that Cat states are highly unstable against entangle- 
ment with environment |9| . Despite such difficulties, Cat 
states on the mesoscopic scale have recently been cre- 
ated 1 1 1 - Since the discovery of BEC, there have been 
suggestions to produce "bigger" Cat states using atomic 
Bose condensates. These include using binary mixtures 
of Bose condensates | [Ll| , as well as performing projec- 
tions on the coherent states [^2|. The attractive Bose 
gas in a double well is mathematically similar to the bi- 
nary mixtures of Bose condensates ]fl| but is a much 
simplier system. As we shall see, it has Cat-like ground 
states over a wide range of parameters, a fact that does 
not seem to be generally appreciated. 

In a separate development, recent studies of spin-1 
Bose gases with antiferromagnetic interaction show that 
strict spin conservation will lead to fragmented ground 
states, which can be Fock like (with zero spin fluctuation) 
or superfragmented (with N 2 spin fluctuation), depend- 
ing on the magnetization of the system || . Although the 
superfragmented states of spin-1 Bose gas and attractive 
Bose gas in double wells assume very different forms, the 
origin of their formation and their fluctuations are iden- 
tical. That it is caused by the existence of degenerate 
minima in the interaction energy in number space and 
the quantum fluctuations between them. In contrast, a 
Fock state is caused by a single deep minimum in the in- 
teraction energy which strongly suppresses number fluc- 
tuations. This mechanism can be seen in spin-1 Bose gas, 
and is best illustrated in the double well example below. 

The Two Site Boson Hubbard Model: Consider 
the Hamiltonian 

H = -t{a)b + tfa) + - [n a (n a - 1) + n b {n b - 1)] , (1) 

with n a + nb = N, where a and b destroy Bosons at site 
(i.e. well) a and b, n a — a^a, — b^b, N is the total 
number of Bosons, t > is the tunneling matrix element, 
and U is the interaction between particles which has the 
same sign as the s-wave scattering length p3| . Eq.(|l|) is 
meant to describe the actual system when reduced to the 
lowest doublet (a ± b) . Although there are residual terms 
in the effective Hamiltonian in the reduction, and that 
t and U in eq.(|l|) have density dependence, we shall not 
consider such features because they do not affect the ba- 
sic physics of the problem (i.e. the competition between 
tunneling and interaction) . Instead, we shall focus on the 
simple model (eq.(|l])) which deserves to be studied in its 
own right. 
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For a system with N Bosons, the Hilbert space is 
{\£) = |f + £,% - £)}, \e\ < N/2, where \N a ,N b )= 
rf N «tf N » | Q j l^N a \N b \ i s t he state with N a and N b Bosons 
at site a and b. We shall take A even, a choice of con- 
venience that has no effects on our results. If \^f) = 
J2i ^i\£) is an eigenstate with energy E, eq.(|l|) implies 

E^ e = - (t t -x^t-i + tt9 t+1 ) + U£ 2 ^ e , (2) 



where t t = Jf (f + 1) - £(£ + 1). It is clear that 
tg strongly favors large amplitudes near £ = 0. This 
is reflected in the non- interacting ground state, |C) = 
(2 N N\)- 1 / 2 {a)+tf) N \0), which shows that if is a Gaus- 
sian centered at I = with a width er c = \/N/2 since 
j^^jVig^/JV for ^ >> L when [/ > o, the 

potential U£ 2 suppresses particle fluctuations between a 
and b, narrowing the Gaussian toward the delta-function 
S e=0 , (corresponding to the Fock state \N/2, N/2)) [0 |§. 
The ground state for general U > (referred to as the 
"repulsive family" ) is 

-OM, h-{h + wT (3) 

where G a {£)= e~ e V 2 " 2 / '(na 2 ) 1 / 4 . EqM follows from the 
fact that in the continuum limit eq.(0) near £ = is the 
Schrodinger equation of a simple harmonic oscillator. Its 
validity can also be verified numerically. 

For U < 0, the potential Ui 2 favors large amplitudes 
at I = ±A/2. It competes strongly with tunneling and 
tends to split the Gaussian tyf at U — into two. These 
features are found in the numerical solutions of eq.((|), 
which also show to a good approximation, the ground 
state (referred to as the "attractive family") is 

¥f ] = Q(*f + ¥?) (4) 

where Q is a normalization constant, ^ R ^= G a (£ — A), 

ty^Li • (See figure 1). In other words, the ground 
state is a superposition of the coherent/Fock-like states 
\L) and \R), |*<-)> = Q(\L) + \R)), 

\L)=J2 G °(<iMn + ,n_, \R) = J2 G -(9)\q)N_,N + , 
9 q 

(5) 

where \q) N+ , N ^_ = \N + + q, N_ - q) , N± = § ± A, and 
N + is the average number of a (or b) particle in \L) 
(or \R)). We also found numerically that the quanti- 
ties A* = A/ (N/2) and a* = a/a c (a c = ^JW/2) are 
functions of the combination UN/ 1 only jl4| , with a di- 
viding behavior at UN/t = — 2 as shown in figure 2. 
As UN/t increases below —2, we have A — > N/2 and 
a — * 0. The system is driven towards the "extreme" Cat 
state |Cat*) = (\N,0) + \0,N))/V2. In fact, the over- 
lap between \L) and \R) vanishes rapidly when A > a, 



or A* /a* > y/2/N. For a system with N = 1000 par- 
ticles, we find that \L) and \R) cease to overlap when 
UN/t < —2.1. The system is essentially a superposition 
to two non overlaping mesoscopic condensates. 

Fragmentation and Super-fragmentation: That 
interaction will lead to fragmentation can be seen from 
that the single particle density matrices of both |1^ + )) 

^i* < - , )./-(J11 ( A ) ) = f('i).^ 

x = e~ 1 /( 4<T ) for I'I'M) with a given in eq.(^), and x — 
(e-VC^ 2 )^ _ (2A)2 +e -(A+i/2)V- 2 ) /(i + e- A2 /° 2 ). for 

l^^ - )) with a given in fig. 2. The eigenvalues of p are 
A± = (N/2)(l ± x). When U -> 0, we have x = 1 for 
both signs of U, and p has only one macroscopic eigen- 
value (A+ = N and A_ =0). This is expected since the 
ground state reduces to the coherent state |C). As \U\ 
increases, x — > for both signs of U. The condensate 
becomes fragmented since both eigenvalues of p becomes 
macroscopic, A + ,A_ — > N/2. 

However, when examining the internal number fluctu- 
ations AN 2 =((N a - (N a )) 2 ), one finds AN 2 = a 2 /2 for 
which vanishes as U/tN increases. In contrast, 

AN a= X+ 1+e -tIw fOT Which iS ° f OTder N2 

since A ~ N and a - y/N for large \U\N/t. The re- 
sults of AiV^ obtained from numerical solution of eq.(||) 
are shown in figure 2. We have plotted ln(4AA r ( ^)/lniV 
instead of AA^ because the former assumes the simple 
value of 2 and 1 for the "extreme" Cat state |Cat*) and 
the coherent state |C) respectively. Fig. 2 shows that 
AJV^ reaches a substantial fraction of its maximum value 
A 2 /4 over the interval A(UN/t)= (-3,-2). Thus, from 
the viewpoint of achieving a superfragmented structure, 
it is not necessary to go deeply into the Cat regime, (i.e. 
\UN/t\ » 1). 

At first sight, such an interval may seem physically ir- 
relevant because for any finite interval A(UN/t), the cor- 
responding range in U/t vanishes as 1/N in the thermo- 
dynamic limit, and that the system is either deep in the 
Cat regime or is a coherent state depending on whether 
U < or U = 0. This, however, is not true. The reason 
is that in general, the tunneling matrix clement t de- 
pends on an energy barrier V in an exponential fashion, 
\nt oc —V a ; (and V is proportional to the intensity of the 
Laser producing the barrier). Thus, the range AV corre- 
sponding to the interval A(UN/t) is only proportional to 
InjV, making the system highly tunable from one regime 
to another. Two other facts also make this transition 
region relevant. Firstly, quantum gases are mesoscopic 
instead of macroscopic systems, with N < 10 5 instead 
of N ~ 10 23 . Secondly, recent experiments have shown 
that the scattering length of Rb 85 can be tuned to zero 
by varying the external magnetic field p5[ . These show 
that it is possible to have a system with a small UN 
even for N as high as 10 5 . The fact that the ground 
state changes continuously from a coherent to a Cat-like 
structure over a wide range of parameter allows us to 



2 



explore how the phase coherence of the system as the 
crossover takes place. 

Interference of Attractive Bose Gas in the Su- 
perfragmented regime: As mentioned earlier, many 
authors have pointed out that there will be no distinc- 
tions between the interference patterns of a coherent 
state and a Fock state. The numerical evidence of this 
effect was given Javanainen and Yoo (JY) ||. Later, 
using analogies with quantum optics, Castin and Dal- 
ibard (CD) simulate the particle collection process by 
the "beam splitter" operators a ± b and showed explic- 
itly how the measurement process changes a Fock state 
into a coherent state. The exact spatial pattern, how- 
ever, was not derived. In the following, we shall modify 
the calculation of CD to obtain the spatial interference 
pattern of the attractive family eq. (Q) . We shall see that 
the operators for particle collection are different from the 
"beam splitter" operators, and that our calculation when 
applied to the Fock state furnishes a derivation of JY's 
result 1 5). 

To illustrate the key features, it is sufficient to con- 
sider the one dimensional case. At time t = 0, the trap 
is turned off and the condensates at a and b begins to 
expand and to overlap. For simplicity, we shall assume 
the atoms expand as non-interacting particles for t > 0, 
which implies i>(x,t) = ds e iM( - x ~ s ? ^ 2ht ip(s). If 



a and b are Wannier states localized at ±x , we then have 
. e -iC{x)/2 b ^ = (2x M/ht)x, and 
,iM{x 2 +x )/2ht _ jf B osons were photons, 



■/ 2 n hi 



7 = ^ 

the operators (a+b) k+ (a—b) k - represent detecting differ- 
ent number of photons in different beam splitters. How- 
ever, in an interference experiment, particle detections 
are products of t[)(x,t), which are specific combinations 
of a and b rather than than products of (a+b) k+ (a — b) k - . 

Next we consider a series of particle detectors located 
at Xi,i = 1, to D. The joint probability of detecting a 
total of k particles (k << N) with ki particles in the 
detector at X4 is, (see Appendix), 



{N-k)\ k\ 



N\ 



(6) 



where {kf } and {kf} are the optimal set of Vl and Vr 
respectively. 
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To calculate {k ■ }, we consider the coherent state 
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■jtf) N \0) 



(J) 



where u = e~ la / 2 cos^, v = e ta ^ 2 sm§, cos 2 f = N+/N. 
It follows from eq.Mthat (a^a) = N + , {tfb) = N-, 
N± = N/2 ± A. Eqjjfy has the expansion 



N 



i(A+q)a\~ 



\Q)n + ,1> 



(8) 



where V g o) = G ao {q) with ct 2 = 2N+N^/N up to 1/N+, 
l/N- corrections. Inverting eq.(|J), we can express 
\q)n ,n_ and hence \L) [eq.(||)] in terms of |a,/3)jv. Wc 
then have 



6|i)=r^/(g) f %^ A +*f[[W Xi (a)] ki \a,P) 

q • ) ^ 7T 77 i=l 



N-k 



(9) 

G a (q)/G a Jq), W x (a)= e^/ 2 u + 
Since \a,j3)jsr-k is close 
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where f(q) 

e -,C(x)/2 U; andr = 7 Y__. 

to the Fock state e lAa \N+, N-), and since the range of 
q is restricted to I/er, the combination e ia ^ A+q '\a, (3) N-k 
varies slowly with a. Since fcj >> I, one can use the 
method of steepest descents to determine the phase an- 
gle a* that maximizes the magnitude Y\a=x I^k«(oj)| % 
and we have 0\L) oc \a*). This shows that the measure- 
ment process (specified by the set {h}) projects the state 
\L) into the coherent state \a*,/3). 

To calculate Vl, we note that (a',(3\a,P) w 
exp[-(^sm 2 /3)(a' - a) 2 +i^{N ~ k)(a' - a)] for 
N-k » 1. We then have 

da 



V L {{h}) = i 1 k 



D 



2tt 



i/»i 2 n 



\W Xi (a)\ 
hi 



2k i 



(10) 



where O — Ylii* \ x i) removes ki particles at n, \^>)n 
is a normalized state with N particles, and Y^i=i ^« = ^- 
The measured density at Xi is given by the most proba- 
ble set {ki} which optimizes V({ki}), i.e. n(xj) — kj. In 
case there are many such sets, the measured density will 
change from experiment to experiment, as each experi- 
ment samples a different optimal set. 

For Cat like states (eq.(0)), we have V({ki})= 
Q 2 [V L ({h}) + V R {{ki})+V LR X{ki})}, where V L and Vr 
are eq.(g) evaluated at |^) = \L) and \R), and Vlr 
is eq.(||) with the norm replaced by {L\0^0\R) + c.c. 
Since Vlr depends on the overlap of \L) and \R), it is 
non- vanishing only within the range of UN/t such that 
A < a, For A > a, V is dominated by Vl and Vr, and 

£ 

the interference pattern n(xj) is proportional to fc_- +fc_- , 



where f(a)= E 9 /(?)e M9 
where <7 ( , 2 ; rr , 2 2 



a , and 77 



As the 



ground state becomes more Cat-like, a — > 0, f(a) has 

a weak a dependence. To find the optimal set {kf}, 
we rewrite the kj. in eq.( |lO|) u sing Stirling formula, and 
optimize the product in eq.(|10|) using method of steepest 
descent subject to the constraint ^ ki = k. One then 



obtains k i oc \W X% (a*)| 2 , or 



kf = A 



I + sin/Jcos 



/ 2Mx a x 

\ ht 



(11) 



where A is a constant. The stationary condition for a* 
can be derived in a straightforward manner and will not 
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be presented here. Repeating the calculation for Vr we 

2j R 

find k i — k i . The measured density n(xj)= 
therefore consists of a uniform background and a sinu- 
soidal oscillation with wavevector hT 1 M(2x a /t). The 
amplitude of oscillation sin/3 vanishes as the systems be- 
comes more Cat-like, since (3 — > as A — » iV/2. This is 
in contrast to the repulsive case, where the interference 
pattern is independent of the barrier height which causes 
the system to fragment into a Fock state. 

Appendix. Eq.(||) was derived in the Appendix in 
H using continuous quantum measurement theory. It 
can also be derived from the following elementary con- 
siderations: For a complete set of states (labelled by p) 
with creation operator {^4p}, the probability of detect- 
ing a particle in state i in an iV-particle system \^f) is 

^ (1) = 4^ = U^lf/mm 2 )- The probabil- 
ity of detecting a second particle in state i after the 
first particle is detected (provided that the state evolves 
very little between detections) is V {2) = WA^A^))]] 2 
/[(7V-1)||A 2 |*)|| 2 ]. The joint probability of detecting k t 

particle in state * is vf ] vf ~ 1] ..v\ 1] = ^#PM*)|| 2 
/||<I')|| 2 . Eq.(^) is obtained by generalizing to the detec- 
tion to more than one states and with Ai replaced by 
if){xi). The combinatoric factors in ki accounts for the 
arbitrary ordering in the detection of different i states. 
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Caption 

Figure 1. Numerical solution of calculated from 

eq.(g) for different UN/t for a system with N = 1000 
particles. The results for U < can be well fitted by the 
functional form in eq. (Q) with parameters A and a shown 

in figure 2. For UN/t < — 2, V E'^~' ) begins to split up into 
two Gaussians. The split-up is complete for UN/t « 
-2.1. 

Figure 2. With the solutions of eq.(Q) for U < well 
fitted by eq.(Q), we find that for different (N,U, and 
t), each of the parameters A* = A/ (N/2), <j/ct c , and 
111(4 AN 2 ) /\nN when plotted against UN/t falls into a 
single curve. Note that ln(4AiV 2 )/lniV= 2 and 1 for the 
Cat state (\N, 0) -I- |0, N))/V2 and the coherent state \C) 
respectively. At UN/t = -2, \n(4AN 2 )/\nN = 1.32. 
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